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ABSTRACT 



We recently discovered a large number of highly active Be stars in the open cluster NGC 
3766, making it an excellent location to study the formation mechanism of Be star disks. To 
explore whether similar disk appearances and/or disappearances are common among the Be 
stars in other open clusters, we present here multiple epochs of Ha spectroscopy for 296 stars in 
eight open clusters. We identify 12 new transient Be stars and confirm 17 additional Be stars 
with relatively stable disks. By comparing the Ha equivalent widths to the photometric y-Ha 
colors, we present a method to estimate the strength of the Ha emission when spectroscopy is not 
available. For a subset of 128 stars in four open clusters, we also use blue optical spectroscopy 
and available Stromgren photometry to measure their projected rotational velocities, effective 
temperatures, and polar surface gravities. We combine our Be star detections from these four 
clusters to investigate physical differences between the transient Be stars, stable Be stars, and 
normal B-type stars with no line emission. Both types of Be stars are faster rotating populations 
than normal B-type stars, and we find no significant physical differences between the transient 
and stable Be stars in our sample. 



Subject headtngs: stars: emission-line, Be - 
3766, NGC 4755, NGC 6231) 



■ open clusters and associations: individual (NGC 3293, NGC 



1. Introduction 

Classical Be stars are a class of B-type stars 
with circumstellar disks that cause Balmer and 
other line emission. As a population, they rotate 
faster than their nonemission, B-type counterparts 
with a rotational velocity comparable to the crit- 
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. There are three pos- 



sible reasons for the rapid rotation of Be stars: 
they may have been born as rapid rotators, spun 
up by binary mass transfer, or spun up during the 
main-sequence (MS) evolution of B stars. The ro- 
tation is likely combined with weaker processes, 
such as nonradial pulsations (NRPs) or magnetic 
fields to move material from the stellar surface into 
the disk ( [Porter fc Rivinius||2003 1 . 

We began studying the Be phenomenon in open 
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clusters with a photometric and spectroscopic in- 
vestigation of NGC 3766 ( [McSwain fc Gies|2005a 
hereafter Paper 1). We then appUed our photo- 
metric technique to perform a survey of 55 south- 
ern open clusters to examine the Be star ages 
and dependence on main-sequence and post-main- 
sequence evolution ( McSwain fc Gies|2005b here- 
after Paper 2). With 52 definite Be star detec- 
tions and 129 Be candidates identified in that sur- 
vey, our results indicated that the fraction of Be 
stars in a cluster increases with age until 100 Myr. 
Be stars are most common among the brightest, 
most massive B-type stars above the zero-age MS 
(ZAMS). A spin-up phase at the terminal-age MS 
(TAMS) cannot produce the observed distribution 
of Be stars, but up to 73% are consistent with 
spin-up by mass transfer in a close binary system. 
Most of the remaining Be stars were likely rapid 
rotators at birth. Other recent surveys of Be stars 
in Galactic and Magellanic Cloud open clusters 



have found similar results ( 


Wisniewski fc Bjork- 


man||2006 Mathew et al.||2008 


). 



During a foUowup spectroscopic study of the 



cluster NGC 3766 ( |McSwain et al.||2008| here- 
after Paper 3), we discovered unprecedented ac- 
tivity among its Be stars. Between 2003-2007, 
we observed disk appearances and/or near disap- 
pearances in 11 of its 16 Be stars. The observed 
disk formation rates implied a very slow equato- 
rial surface flow on average, but the surface flow 
can be an order of magnitude faster at times. Our 
measurements of the changing disk sizes in NGC 
3766 are consistent with the idea that transitory, 
nonradial pulsations contribute to the formation 
of these highly variable disks. 

One of the primary motivations of this work 
was to confirm many of our Be candidates from 
Paper 2 with spectroscopic observations and im- 
prove the fraction of known Be stars in southern 
open clusters. We also sought to discover addi- 
tional Be stars that exhibit frequent changes in 
their disk states, transitioning between the Be and 
"normal" B-type spectra, as we observed among 
most of the Be stars in NGC 3766. Finally, spec- 
troscopic observations allow us to probe the ro- 
tational velocities and other physical parameters 
that may distinguish Be stars from non-emitting 
B-type stars. 

Throughout this work, we adopt the term 
"transient Be stars" to describe any Be star ob- 



served in both the Be and normal B-type spectro- 
scopic states. Likewise, we use the term "stable Be 
stars" to describe Be stars whose disks are present 
in all observations, even though the strength of the 
disk emission may change dramatically. 

Here we present multiple epochs of Ha spec- 
troscopy of eight open clusters, and we identify 12 
previously unknown transient Be stars from this 
sample. We also confirm 17 additional Be stars 
that have relatively stable disks. For a subset 
of four open clusters, we also present blue opti- 
cal spectroscopy and measurements of the effec- 
tive temperature, T^s, surface gravity, log g, and 
projected rotational velocity, Vsini. Combining 
these results with our previous measurements of 
members of NGC 3766, we compare the total sam- 
ples of transient Be stars, stable Be stars, and nor- 
mal B-type stars to explore physical differences in 
these three populations. Other than the tempo- 
rary disappearance of their disk emission, we find 
no physical differences between transient Be stars 
and the Be stars with relatively stable disks. 

2. Observations 

We obtained red optical spectra of 296 Be and 
B-type stars in 8 open clusters during observing 
runs in 2005 February, 2006 May, and 2007 May 
using the CTIO Blanco 4-m telescope with the Hy- 
dra multi-fiber spectrograph. For a subset of 4 
open clusters, we also obtained blue optical spec- 
tra during the 2006 and 2007 runs. The spec- 
trograph setups differed slightly between all these 
runs, and the details are summarized in Table [T| 
We used slit plates and/or order-sorting filters in 
some of our observing runs, and when either was 
not used. Table [T] contains a null value. 

We selected the targets for each run by giving 
highest priority to the known Be stars in the clus- 
ter. We then selected other B-type stars in the 
cluster by ranking them according to their y-Ha 
color to preferentially select any weak emission. Be 
star candidates from our photometry (Paper 2). 
In Paper 3, we published results from 47 mem- 
bers of NGC 3766, and we here we expand that 
sample with red and blue spectra of 10 additional 
members observed in 2007 May. Also in Paper 3, 
four Be stars did not have available blue spectra 
to measure their Vsini, and we included them in 
the new target list. For the other clusters, we used 
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the same fiber assignments to achieve a target list 
as consistent as possible, but a broken fiber pro- 
hibited a few stars from being observed in later 
runs. Also, the 2007 blue spectra suffered from 
much lower signal-to-noise ratios (S/N) due to the 
slit plate inserted in front of the fibers (possibly 
misaligned and blocking light), so several targets 
from this configuration were excluded from our 
measured sample. 

We generally began the Hydra observations by 
taking short exposures and then parking the fibers 
used for the brightest stars to avoid saturation in 
longer exposures of 1800 s. Therefore we usually 
obtained up to 5-7 exposures of each star. We also 
observed a HeNeAr comparison lamp source just 
before and after each set of cluster observations 
for wavelength calibrations. 

All of the spectra were zero corrected using 
standard routines in IRAF, and they were flat 
fielded, wavelength calibrated, and sky subtracted 
in IRAF using the dohydra routine. Of the 7 clus- 
ters observed, only NGC 3293 and IC 2581 had 
widespread nebulosity that was difficult to remove 
completely during the sky subtraction. The effect 
was pronounced in the Ha line in our red spectra, 
but the nebulosity presented much less of a prob- 
lem in the H7 and other Balmer lines in our blue 
spectra. For each set of Hydra spectra, we trans- 
formed the observations to a common heliocentric 
wavelength grid, co- added them to achieve good 
S/N for each star, and rectified the total spectrum 
to a unit continuum. 

3. Be Star Variability and Ha Equivalent 
Widths 

In Paper 3, we reported the discovery of 11 
transient Be stars in NGC 3766 whose disks ap- 
peared and/or nearly completely disappeared dur- 
ing our 4 years of observations. Only five of the 
observed Be stars had relatively consistent disks 
during the same period. It turns out that NGC 
3766 is not the only cluster containing such vari- 
able Be stars; during two or three epochs of red 
spectroscopy, we have discovered many new tran- 
sient Be stars in nearly every cluster reported here. 
The Ha profiles of each Be star are shown in Fig- 
ures [THSl 

We measured the equivalent widths of the Ha 
line profiles, Whu, by directly integrating across 



each rectified spectrum over the range 6535-6591 
A to include the broad wings expected for this 
strong line. The error in Wna is about 10% due 
to noise in the continuum regions. We use the 
convention that a pure absorption line has a posi- 
tive equivalent width. Many authors prefer to in- 
vestigate Be emission strength by subtracting the 
photospheric absorption line profile and measur- 
ing Whoi from the net emission profile only. We 
avoid that here for two reasons: we would like to 
compare the raw Wna for both emission line and 
non-emission line stars to their photometric colors, 
and we do not know an accurate photospheric line 
strength without knowing the effective tempera- 
ture and gravity of the stars. We measure these 
stellar properties below for a subset of our targets, 
but for consistency, we present only the raw Wna 
for all stars. 

The resulting Wna are listed in Table |2] 
columns 2-4. In column 5, we also identify the 
Be stars and candidates found in our photometry 
(Paper 2). Finally in column 6, we list the Be star 
status based upon our red spectroscopic observa- 
tions only. While the photometric method is often 
unsuccessful identifying small or transient disks, 
it is usually very accurate when identifying strong 
emitters. Only three stars were identified as Be 
in Paper 2 that are not confirmed here: Collinder 
272-565, NGC 4755-43, NGC 6231-67. We argue 
below that the first is likely a foreground star; the 
latter two may be Be transients that have not yet 
been confirmed spectroscopically. 

Several Be stars exhibit dramatic changes in 
their disk states, alternately showing strong disk 
emission and a near complete disappearance of the 
disk: CoUinder 272-562, Collinder 272-621, Hogg 
16-185, NGC 3293-83, NGC 4755-18, and NGC 
4755-31. Normal B-type stars are not expected 
to have strong stellar winds with emission form- 
ing in the Ha line, so any weak emission partially 
filling in the line or the wings is more likely due 
to a very small Be disk. Several stars show line 
profile variations consistent with weak disks dur- 
ing one or more observations: Collinder 272-767, 
NGC 3293-74, NGC 6231-2, NGC 6231-90, and 
NGC 6231-221. Finally, the star IC 2581-113 ap- 
pears to have a transient shell profile, indicating 
absorption through an edge-on disk. We classify 
all 12 of these stars as Be transients. Seven of 
these transient Be stars have available blue spec- 



3 



tra, and combined with the 11 transient Be stars 
observed in NGC 3766, we have a large sample to 
compare the transient population with stable Be 
stars, discussed below. 

An additional 17 Be stars exhibit line emission 
from their disks that are present in all of our ob- 
servations, yet these "stable" disks vary in Who 
just as much or more than the transient disks. In 
Paper 3, we found very similar mass gain and mass 
loss rates for the disks among all of the Be stars, 
suggesting similar origins for all of the disks. The 
Be stars presented here have variability properties 
that are highly consistent with those observed in 
NGC 3766. 

In total, we have identified 29 Be stars from 
our Ha spectroscopy of CoUinder 272/Hogg 16, 
IC 2581, NGC 3293, NGC 4755, NGC 6231, and 
NGC 6664, as well as 16 Be stars in NGC 3766 
(Paper 3). Based on multiple epochs of Ha spec- 
troscopy for most of these 45 Be stars, we present 
lower limits of each cluster's transient Be star frac- 
tion in Table |3] Here, we calculate the percentage, 
p, of transient Be stars, t, compared to the total 
number of Be stars, n, using 



P - 



t 

n 



(1) 



Although the total populations of Be stars in each 
cluster are relatively small, the standard error as- 
sociated with this fraction of transients is approx- 
imately given by the root of the variance for a 
binomial distribution. 



p{l -p) 



(2) 



For NGC 6664, this standard error should be zero 
since p = 0. However, we find a more realistic 
error by allowing p to vary until dp ~ p. The 
Poisson errors associated with t and n are ^/i and 
i/n respectively. The clusters CoUinder 272/Hogg 
16, NGC 3766, and NGC 6231 have the highest 
rates of variable Be star disks. 



Meilland et al. ( 2006 1 discuss two mechanisms 



that may account for our observed transitions from 
the Be to B spectroscopic phases: discrete mass 
loss outbursts that form expanding ring-like struc- 
tures around the star, and continuous mass loss 
that slowly decreases in magnitude and allows the 
disk to dissipate. Based on our measurements of 



the changing disk masses in both stable and tran- 
sient Be stars (Paper 3), we favor the first sce- 
nario. In that work, we measured mass gain rates 
for the disks that appear to increase by an order of 
magnitude during short time intervals. However, 
our spectra are not of sufficient resolving power 
to prove whether the emission lines evolve accord- 
ing to the ring scenario. Higher resolution spectra 
will determine whether the peak separation in the 
emission lines remains constant while the high ve- 



locity wings disappear, as Meilland et al. (20061 



predict for such a discrete mass loss scenario. 

In Paper 1, we showed that in NGC 3766, the 
observed W^a ^^nd y—Ha colors have an approxi- 
mately logarithmic relationship according to 

y — Ha = C 2.51og , (3) 

where FWHM is the range covered by the nar- 
rowband Ha filter and C is a constant. Using our 
mean observed Whoi presented in Table [2] and the 
y— Ha colors from Paper 2, we have performed a 
similar fit for each cluster presented here. The 
resulting fitting constants C for each cluster are 
summarized in Table |4] and the adjusted colors 
y— Ha — C for all stars with measured Wna are 
plotted in Figure |4j The observed ranges in Wru 
are shown with horizontal error bars. Since Wua 
is highly variable even among the stable Be stars, 
most of the scatter in this plot is explained by 
the time separation between our photometric and 
spectroscopic observations. Some additional scat- 
ter is expected since the disks likely vary in den- 
sity, opacity, and inclination. The outlying star 
with WHa = 10 A (CoUinder 272-565) is either a 
transient Be star or a foreground star (more likely 
since transient Be stars usually have bluer y— Ha 
colors due to their smaller, temporary disks). 

Since reddening affects the continuum flux mea- 
sured over the intermediate-band y fllter more 
than the narrow-band Ha filter, the constant C 
is dependent on the cluster's reddening, E{b~ y). 
It will also depend upon the absolute scale of the 
Ha magnitude. We show in Figure [5] that for our 
calibration of y— Ha, there is a direct relationship 
between C and E{b — y). We find a linear fit 

C = (0.97 ± 0.12)£;(6 -y) + (0.14 ± 0.04) (4) 

between the two parameters. Since spectroscopy is 
often unavailable in many observing programs, it 
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is thus possible to estimate Whu of a Be star from 
its photometric properties y—Ha and E(b — y). 
A simple radiative transfer model such as that of 



Grundstrom & Gies (2006) may then be applied 
to investigate the physical size of the disk. How- 
ever, due to the large scatter in the observed rela- 
tions we caution that absolute disk measurements 
are very imprecise and only differential changes 
in ?/— Ha should be applied to estimate first-order 
changes in the disk's extent. Since the Ha filter 
system is not universally standard, different filters 
or photometric calibrations may produce a slightly 
different relationship between C and E{b — y) that 
must be determined before estimating Wna- 

4. Physical Parameters of Stars 

In Paper 3, we presented the physical param- 
eters of 38 Be and B-type stars in NGC 3766 
measured from blue optical spectra. Here, we 
present measurements for additional members of 
NGC 3766 and three other clusters, for a total 
blue sample of 128 Be and B-type stars. The blue 
sample includes only 17 of the transient Be stars 
and 12 of the stable Be stars, and we discuss the 
measurements of V smi, T^s, and log i^poiar below. 

4.1. V smi Measurements 

Our technique for measuring V sin i for B-type 
stars is described in detail in Paper 3. In short, we 
generated a grid of model spectra using SYNSPEC 



( |Lanz fc Hubeny||2003[) and Kur ucz ATLAS9 at- 
mospheric models ( Kurucz||1994 1 assuming plane- 
parallel geometry and local thermodynamic equi- 
librium (LTE). We also adopted solar abundances 
and a microturbulent velocity of 2 km s~^ for all 
stars. For each star, we made a preliminary es- 
timate of its effective temperature and gravity. 
Toff and log g respectively, by comparing the ob- 
served H7, H5, and He I line profiles to our grid 
of Kurucz spectral models. To measure V^sin?, 
we compared the observed He I A4388 and He I 
A4471-|-Mg II A4481 hue profiles to the model pro- 
files convolved with a limb-darkened, rotational 
broadening function and a Gaussian instrumen- 
tal broadening function. We determined the best 
fit over a grid of values by minimizing the mean 
square of the deviations, rms^. The formal error. 
Ay sin i, is the offset from the best-fit value that 
increases the rms^ by 2.7rms^/iV, where N is the 



number of wavelength points within the fit region. 
Our measured V sin i and Ay sin i are listed in 
columns 3-4 of Table O 

In nearly all of the Be stars in our sample, the 
He I lines do not exhibit obvious signs of emission. 
These lines arc, in general, much less susceptible 
to emission than the H7 or H(5 lines. Nevertheless, 
we consider our V sin i measurements for Be stars 
to be lower limits since the He I lines may contain 
some weak emission in Be stars, partially filling 
and narrowing their line profiles. On the other 
hand, those clusters with strong nebular emission 
may have weak He I emission that may cause these 
lines to appear more shallow, hence resulting in an 
overestimate of V smi. 

We compared our measured V sin i to the 
measurements available for each cluster in the 
WEBDA databastj5 The agreement is gener- 
ally good, although we find a large scatter due 
to the systematic differences among authors. For 
18 stars, WEBDA lists V smi measurements from 
multiple sources. The differences between our 
measurements and WEBDA is comparable to the 
differences between the various sources. We mea- 
sured the standard deviation, cry sin i, of the dif- 
ference between our measurements and the liter- 
ature to be 25.7 km s~^. Stars with more than 
3o'y sin i difference in V sin i are probably unre- 
solved double lined spectroscopic binaries (SB 2s), 
and we find three new SB2 candidates in our sam- 
ple: NGC 3293-27, NGC 6231-132, and NGC 
6231-137. We note that IGarcia fc MermilliodI 



(|2001j) argue that NGC 6231-132 is an SBl with 
a period of 8.9 d and that NGC 6231-137 is a /3 
Cep star and a probable SBl. These and other 
spectroscopic binaries identified by blended lines 
in our spectra or by classification in the literature 
are identified in column 15 of Table |5] 

4.2. Tcff and logy Measurements of B stars 

We measured Tcff and log g for each B-type 
star using the technique described in Paper 3, 
and we summarize the procedure here. For the 
coolest stars with Toff < 15,000 K, we used the 



"virtual star" method of Huang & Gies ( 2006 1 , 



which assumes a spherically symmetric geometry 



^The WEBDA database is maintained by J.- 
C. Mermilliod and is available online at 
Ihttp: //obswww.unige.ch/webda/navigation.htmll 



and constant Tes and logg across the stellar sur- 
face. They generated detailed H7 hne profiles 
using line-blanketed, LTE Kurucz ATLAS9 and 
SYNSPEC codes. Huang & Gies show that the 
H7 wing strength and line equivalent width can 
be used as starting parameters in a line profile fit 
to obtain unique values of Tcft, logg, and their 
corresponding errors. 

To measure the hotter B-type stars, we used 
the new TLUSTY BSTAR2006 grid of metal line- 
blanketed, non-LTE, plane-parallel, hydrostatic 
model spectra ( |Lanz fc Hubeny 20071. We used 
their models with solar metallicity and helium 
abundance and a microturbulent velocity of 2 
km s~^. We compared the H7 line profile to the 
rotationally and instrumcntally broadened model 
spectral line profiles at each value in the grid, min- 
imizing rms^ across the line region. We then re- 
fined our measurements to a higher precision using 
a linear interpolation between the available line 
profiles in the grid. Finally, we determined the 
errors, ATcff and A log 5, from the values which 
produce a rms^ no more than 2.1 rms^ / N greater 
than the minimum rms^. Our measurements of 
Tcff and log g, with their corresponding errors, are 
listed in columns 5-8 of Table [s] 

For stars with Teff near 15,000 K, we measured 
the parameters using both methods. Since the 
method of Huang & Gies can be applied to the 
hotter stars (although it does not include non- 
LTE effects), it is an appropriate tool to compare 
with TLUSTY in this temperature regime. We 
find very good agreement from both methods, and 
we adopt the TLUSTY results here. 

A few stars in NGC 3293 were contaminated by 
nebulosity in the H7 line, manifested as a narrow 
emission profile in the line core. To measure Teg 
and log g for these stars, we omitted the center 5- 
6 A of the line profile in our fits and used only the 
broader wings. 

Rapidly rotating B stars may be distorted into 
an oblate spheroidal shape, so the surface grav- 
ity at the equator can be much lower than at the 
poles, and this \og gpoiar is a better indicator of 
the evolutionary state of the star. |Huang fc Gies| 
(20061 performed detailed spectroscopic modehng 



of such distorted rotating stars to determine a sta- 
tistical correction factor for log g, averaged over 
all possible i, for a variety of stellar models. We 
made a bilinear interpolation between their mod- 



els to convert our measured log g to log gpoiar for a 
more accurate comparison between slow and rapid 
rotators. 

For each B star, we also measured its mass. 
Mi,, and radius, R-,,, by interpolating between the 
evolutionary tracks for non-rotating stars from 
Schaller et all (|1992|). The errors AM^ and Ai?* 



correspond to our measured ATpff and A log g. For 
simplicity, we assume the polar radius of the star, 
i?p is equal to i?^, and a rotationally distorted star 
has an equatorial radius = 1.5i?p in the Roche 
approximation. The resulting critical velocity is 



(5) 



Our results for logt/poiar, Af*, AAf*, i?^, and Ai?^, 
and Vcrit are also listed in Table |5] columns 9-14. 

4.3. Tcff and log (7 Measurements of Be 
stars 

We measured Toff and log g for the Be stars us- 
ing available Stromgren mi, Cf, and (3 indices for 
our targets, available from the WEBDA database, 
as described in Paper 3. To determine the dered- 
dened indices, we adopted the reddening, E(b—y), 
listed in Table HI We used the relation 



E{b-y) = 0.745 E{B-V) 



(6) 



to relate E{b — y) and E{B — V) when the lat- 
ter was given in the literature. In Paper 3, we 
found that the temperature relation based on the 



dereddened Stromgren indices from Balona (19841 



provides a fair agreement to our measured Teff 
from the TLUSTY model fits for B stars with 
Toff > 15, 000 K. A small correction factor im- 
proves the agreement. We formed a sample of 
temperature calibrators using the B stars in NGC 
3293, NGC 4755, and NGC 6231 with Teff > 
15,000 K (measured from the TLUSTY models 
as described above) and available Stromgren pho- 
tometry. We also included eight calibrators with 
well-known Teff > 15000 K (measured from their 
absolute integrated stellar fiux) from Napiwotzki 



et al. (19931. We performed a linear fit to Teff 
and TBaiona and found a very similar relationship 
to that found in Paper 3. The photometry in 
WEBDA is from a number of different sources that 
often do not agree well, so the scatter in our tem- 
perature relation for these calibrators introduces 
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an error of 1300 K, much higher than we found for 
NGC 3766 in Paper 3. 

To apply this temperature relation to the Be 
stars, we also need an accurate measurement of 
their (3 magnitudes, which are unreliable due to 
their emission at this line. As in Paper 3, we used 
the B star calibrators to find a linear fit to the ob- 
served (co,/3) relation. The mean scatter between 
/3fit and the measured /3 is 0.017, which implies an 
additional error of 424 K in Toff using the corrected 
Teaiona relation above. 

Finally, we measured Toff for all Be stars in 
NGC 3293, NGC 4755, and NGC 6231 with 
Stromgren photometry available in the WEBDA 
database using the adopted /3fit and the corrected 
Tbaiona as described above. Adding the tempera- 
ture errors from both fits in quadrature, we adopt 
a total error of 1400 K. The results for the Be stars 
are listed in columns 5-6 of Table|5] We found that 
all but one of the Be stars have T^s > 15,000 K, 
and the single exception has T^s — 14, 704 K. 
Therefore we did not need to extend the calibra- 
tion sample to include stars with lower tempera- 
tures. 

The values of log g are more strongly dependent 
on /3, so we were reluctant to use our /3fit with 
the Stromgren relation for log g given by |Balona 



(1984 1 . Instead, we used the calculated Teg to de- 



termine the Be stars' bolometric corrections, BC, 



from Lanz & Hubeny (20071, at first assuming 
log 5 = 4.0. We calculated each stellar radius, 
i?^,, and luminosity, L^,, using the measured T^s, 
BC, V magnitude, distance modulus {V — My)q, 
and E{B - V). The adopted {V - Mv)o are listed 
in Table |4j We measured the stellar mass, M^,, 
from the computed Tgs and L^, by interpolating 



between the evolutionary tracks of Schaller et al. 



( 1992 ) . Finally, we obtained a preliminary value 



of logg from and i?^. Since the BC is weakly 
dependent on log 5, we improved the BC from the 
initial estimate and iterated to compute the final 
logg. We adopt a formal error in logg computed 
from ATeff and the quoted errors in {V — Mv)o 
and E{B-V). 

We note that the evolutionary tracks of |Schaller| 
( 1992 1 do not account for rotation, espe- 



et al. 



cially the very fast rotation observed among the 
Be stars. Evolutionary tracks that include fast ro- 
tation reveal differences (sometimes very substan- 
tial) in the luminosity and MS lifetime of rapidly 



rotating stars (e.g., Ekstrom et al. 20081. How- 
ever, the photometric magnitudes and colors of Be 
stars in our work are likely contaminated by en- 
hanced bound-free and free-free emission from the 



circumstellar disks. Carciofi et al. (20061 showed 



that the V magnitude of the Be star 6 Sco in- 
creased by 0.7 magnitudes during a recent disk 
outburst. However, they also found that the emis- 
sion line strength and visual brightness in S Sco are 
sometimes anticorrelated, possibly due to changes 
in the mass loss rate or disk geometry. The uncer- 
tainty in V combined with the large errors in our 
measured Tes for the Be stars in this work negate 
any advantage in using evolutionary tracks for fast 
rotating stars. Therefore we apply the Schaller et 
al. models throughout this work. 

As we found in Paper 3, the agreement is very 
good between the Stromgren logg and log^poiar 
for the nonemission B-type calibrators, so we 
do not perform any further correction to obtain 
loggpoiar for the Be stars. However, several ef- 
fects may cause the Be stars to appear slightly 
more evolved and artificially brightened. Scat- 
tered light and continuum emission from the disk 
can increase the apparent brightness, an effect 
that is nearly impossible to correct without de- 
tailed knowledge of the disk geometry and orien- 
tation. Rotational mixing in these rapid rotators 
can enrich the hydrogen-burning core and increase 



their luminosity during their evolution (Heger & 



Langer 2000| [Meynet fc Maeder|[2000l [Ekstrom et' 
|al.||2008 1. On the other hand, the rapid rotation 
will also decrease the surface temperature at the 



stellar equator (the Von Zeipel effect; von Zeipel 



1924), and depending on the orientation, the fiux 
may decrease. We consider the values of log^poiar 
to be lower limits, and our measurements of R^, 
upper limits, for the Be stars. The final parame- 
ters for these Be stars are listed in Table [H 

5. The Be and B Star Populations 
5.1. Rotational Velocity Distributions 

The cumulative Vsini distribution for all Be 
and normal B-type stars in NGC 3293, NGC 3766, 
NGC 4755, and NGC 6231 are shown in Figure |6] 
and we find that all four are very similar. To test 
the null hypothesis that the distributions differ, 
we used the two-sided Kolmogorov-Smirnov (K-S) 
statistical test with all six permutations of pairs of 
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cluster data. The K-S test indicates that no sin- 
gle cluster has a significantly diff'erent distribution 
of Tasini. The known and candidate SB systems 
may have overestimated measurements of Vsini 
due to line blending, but they do not infiuence the 
relative distributions of the clusters' Vs'mi. We 
measure a mean V sin i = 184 km s~^ with a stan- 
dard deviation of 90 km s~^ for the normal B-type 
stars, including the SB systems. 

A number of environmental factors may influ- 
ence rotational velocities among cluster members: 
the ages, binary fraction, and magnetic fields. 



Huang & Gies (2006 20081 find that the most 



massive B-type stars spin down with age as pre- 
dicted by models due to their radial expansion and 
angular momentum lost by stellar winds. Since 
the four clusters we observed have nearly identi- 
cal ages (Paper 2), we are not surprised that their 
Vsini distributions are very similar. Binary sys- 
tems are expected to spin down faster than single 
star systems due to tidal synchronization processes 
( |Abt et al.|2002| , and [Huang fc Gies|p056l ) found 
a sharp decline in V sin i for binaries in which the 
more massive component has log^poiar 1^ 3.9. The 
cluster NGC 6231 has a very high binary fraction 



of at least 52% among its B-type stars (Raboud 



1996 ) and at least 63% among 0-type stars (Sana 
et al.|2068 1 , so many of our targets are likely bina- 



ries that will lose angular momentum due to tidal 
interactions as they evolve. However, most of our 
targets have log^poiar > 4.0 and do not appear to 
have reached synchronization. Finally, Huang & 



Gies (2006) found that intermediate and low mass 



B stars spin down faster than more massive B-type 
stars, possibly due to additional angular momen- 
tum losses such as magnetic fields. Indeed, two 
magnetic B stars were recently detected in NGC 



3766 ( McSwain 2008 1 . The similarity in the V sin i 



distributions in NGC 3293, NGC 3766, NGC 4755, 
and NGC 6231 suggests that the environments are 
very similar or that the clusters are too young for 
the various environmental factors to produce ob- 
servable differences in rotation. 

Since we find no systematic differences among 
the cluster Vs'mi distributions, we can explore 
trends in V sin i for the combined Be star pop- 
ulations. We have measured Vsini for 12 Be 
stars with relatively steady-state disks (emission 
present in each of our observations) and 17 Be 
stars with significant changes in their disk emis- 



sion state (either outbursts or disk fading events). 
The mean V sin i and its standard deviation for 
the stable (transient) Be stars is 263 ± 66 km s~^ 
(253 ± 74 km s~^). Figure [t] shows the cumulative 
distribution of V sin i for the Be stars with stable 
disks and transient disks compared to the normal 
B-type cluster population. The K-S test indicates 
a relatively high probability (46%) that the tran- 
sient and stable Be stars have the same V sin i dis- 
tribution, so the differences in their distributions 
are not significant. For the combined Be star sam- 
ple, the K-S test indicates only a 0.03% probability 
that their Vsini are drawn from the same distri- 
bution as the normal B-type stars. Both the sta- 
ble and transient Be stars are clearly populations 
of rapid rotators compared to the normal B-type 
stars in the sample, in excellent agreement with 



several other recent studies of Be rotation ( Hunter 



et al.||2008[|Martayan et al.|[2006l [Martayan et aL 



2007 



Paper 3). 



5.2. Tcff and log 5 Distributions 



In Figures [SffTT] we plot the Tog and log g dis- 
tributions for each cluster with the correspond- 
ing evolutionary tracks for nonrotating stars from 



SchaUer et al. ( 1992 ) 



The large population of Be transients with 
15000 < Teff < 20000 K in NGC 3766 is strik- 
ing, but we doubt that this is significant. Due 
to the cluster's older age, its members are more 
evolved and are generally not found at higher tem- 
peratures. In addition, nearly every Be transient 
in the other three clusters lies outside this range, 
and the transients are well distributed among the 
entire observed range of temperatures. The rela- 
tionship between transients and cluster age is not 
clear, and it is worth investigating the transient 
Be star phenomenon in more clusters to determine 
whether there is an evolutionary dependence. 

Like the transient Be stars, the stable Be stars 
are equally distributed over the entire temperature 
range. There is no clear dependence in T^s and 
logg that distinguishes the stable Be stars from 
the transients or normal B-type star populations. 

6. Summary 

From Ha spectroscopy of 296 stars in 8 open 
clusters, we have identified a total of 23 Be stars 
that experience transitions from the Be to B spec- 
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troscopic phase (including those in NGC 3766 that 
we reported in Paper 3). We also find 22 Be stars 
that have relatively stable disks that are consis- 
tently present in our multiple epochs of Ha spec- 
troscopy. Transient Be stars do not appear equally 
common in the open clusters, and we find higher 
rates of disk appearance and disappearance events 
in CoUinder 272, Hogg 16, NGC 3766, and NGC 
6231. These clusters will be important sites to 
study the disk dissipation and formation mecha- 
nisms for Be stars. 

Prom blue optical spectroscopy and available 
Stromgren photometry of 128 Be and B-type stars 
in 4 open clusters, we have measured the T^s, 
loggpoiar, and Vsini of these stars. Our sam- 
ple of physical parameters includes 17 transient 
Be stars and 12 stable Be stars, and we compare 
their rotational properties to find that all of the Be 
stars have a similar distribution of V sin i and are 
more rapidly rotating than normal B-type stars. 
Other than the temporary disappearance of their 
disk emission, the transient Be stars do not exhibit 
significant physical difi^erences from Be stars with 
relatively stable disks. 

We are grateful to the anonymous referee for 
comments that greatly improved this manuscript. 
We thank the staff at CTIO, especially Knut Olsen 
and all of the night assistants, for their hard 
work and patient support during our many ob- 
serving runs. Support for MVM has been pro- 
vided by an NSF Astronomy and Astrophysics 
Postdoctoral Fellowship and by an institutional 
grant from Lehigh University. This material is 
based upon work supported by the National Sci- 
ence Foundation under Grants No. AST-0606861, 
AST-0507219, and AST-0401460. 
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Table 2 
Ha Equivalent Widths 



Cluster, 


HJD (t - 2,450,000), 


HJD (t - 2,450,000), 


HJD (t - 2,450,000), 


Emission status Emission status 


MG ID 


WHa (A) (2005 Feb) 


WHa (A) (2006 May) 


Who, (A) (2007 May) 


in Paper 2 in this work 


CoUinder 272 


3404.832 


3870.178 


4225.757 




42 


8.60 


8.87 


7.35 




44 


7.60 


7.64 


10.22 




157 


8.84 


9.17 


10.45 




188 


8.20 


8.94 


9.11 




199 


10.45 


10.77 


10.94 




222 


10.52 


10.95 


10.97 




225 


9.82 


10.48 


9.29 




254 


4.43 


4.62 


4.63 




267 


5.92 


6.28 


6.23 




288 




10.54 


10.13 




302 


6.54 


6.91 


6.77 




309 


6.62 


7.31 


7.41 




310 


10.66 


10.70 


10.85 




311 


7.74 


7.74 


7.84 




329 


8.37 


9.64 


10.00 




349 


7.99 


8.65 


8.47 




376 


10.79 


10.54 


10.67 




377 


7.50 


8.40 


8.27 




397 


8.93 


9.20 


8.58 




406 




13.00 


12.36 




448 


10.47 


8.92 


9.21 




484 


12.58 


9.13 






508 


-42.93 


-38.08 




Be Be 


513 


-6.90 


-2.73 


-0.46 


Be? Be 


560 


10.33 


10.08 


10.32 


Be? 


561 


9.00 


8.65 


8.62 




562 


-4.56 


-1.70 


2.92 


Be? Be transient 


565 


10.23 


10.52 




Be 


621 


1.01 


4.90 


8.43 


Be transient 


665 


7.08 


7.34 


7.84 




678 




6.79 


7.27 




708 




10.59 


10.45 




717 


6.80 


7.74 


7.62 




755 


11.67 


11.49 


11.44 




767 


7.71 


10.94 


7.52 


Be Be transient 


793 




11.71 


10.97 




Hogg 16 


3404.832 


3870.178 


4225.757 




94 


8.10 


8.26 






125 


4.44 


4.55 


3.93 




134 


12.03 


11.61 






154 


8.47 


8.55 






185 


1.69 


5.54 


2.69 


• • • Be transient 


209 


4.32 


4.54 


4.21 




217 


12.16 


12.14 


14.34 




227 




5.40 






254 




13.17 






317 


4.43 


5.05 


4.77 




ic 2:)«i" 


:!lUi.7U« 




4225.51!)! 




5 


8.72 




11.07 




9 


9.22 




10.04 




10 


5.87 




6.05 


Be? 


19 


6.91 




7.27 




27 


11.12 




11.77 




42 


7.54 




7.21 




44 


5.84 




6.40 
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Table 2 — Continued 



Cluster, 


HJD (t - 2,450,000), 


HJD (t - 2,450,000), 


HJD (t - 2,450,000), 


Emission status Emission status 


MG ID 


WHa (A) (2005 Feb) 


Who, (A) (2006 May) 


Who, (A) (2007 May) 


in Paper 2 in this work 


48 


7.05 




7.21 




73 


12.96 




11.93 




84 
99 


9.46 
10.08 




10.35 
11.34 




105 


10.31 




10.76 




113 


8.40 




10.48 


• • • Be transient 


119 


4.05 




4.34 




122 


4.06 




5.65 




123 


5.80 




7.10 


Be? 


153 


4.53 




4.92 




157 


6.52 




7.11 


Be? 


158 


-14.88 






Be Be 


173 


8.59 




8.43 




179 


7.24 




7.74 




181 


7.09 




7.30 




188 


6.99 




7.83 




210 


2.80 




3.35 




216 


9.10 




9.38 




219 


3.06 




6.90 


Be? 


233 


7.74 




8.57 




243 


—48.51 




—46.44 


Be Be 


247 


7.07 




7.46 




251 


3.69 




3.97 




267 


0.88 




1.05 


Be? 


270 


5.13 




6.64 


Be? 


277 


7.54 




7.55 




295 


3.90 




4.07 




298 






7.21 




314 


8.97 




9.29 




316 


9.32 




9.24 




323 


7.25 




6.27 


Be? 


329 


6.24 




6.34 




370 


9.41 




9.57 




NGC 3293" 


3404.772 


3870.506 


4223.321 




1 


6.02 


5.43 


4.92 




2 


5.47 


5.22 


4.68 




4 


4.24 


3.29 


-0.60 


Be? 


7 


8.10 


7.87 


7.05 




8 


4.75 


4.67 


4.07 


Be? 


9 


5.64 


6.27 






16 


9.08 


9.37 


9.05 




18 


9.27 


8.85 






22 


4.35 


4.43 


4.50 




23 


4.18 


6.14 


6.12 




27 


6.02 


6.26 


6.10 


Be? 


28 


3.96 


4.20 


4.20 




40 


3.51 


4.10 


4.16 




42 


-1.60 


-1.48 


-2.26 


Be Be 


46 


2.14 


2.83 


2.72 


Be Be 


47 


3.32 


3.70 


3.86 




51 


4.58 


5.68 


4.67 


Bo? 


54 


6.57 


5.06 


4.23 


Be? 


61 


6.03 


5.94 


5.27 




62 


3.52 


3.73 


3.88 




67 


6.18 


6.39 


6.12 




74 


2.56 


3.83 


4.16 


• • • Be transient 


83 


0.47 


4.14 


4.15 


• • • Be transient 
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Table 2 — Continued 



Cluster, 


HJD (t - 2,450,000), 


HJD (t - 2,450,000), 


HJD (t - 2,450,000), 


Emission status Emission status 


MG ID 


WHa (A) (2005 Feb) 


Who, (A) (2006 May) 


Who, (A) (2007 May) 


in Paper 2 in this work 


84 


7.22 


7.31 


7.34 


Be? 


90 


4.09 


4.08 


4.27 




95 


3.20 


3.84 


3.97 




99 


— 11.99 


— 15.89 


—20.47 


Be? Be 


100 


3.35 


3.26 


3.65 




107 


7.74 


8.78 


9.01 




110 


3.89 


5.19 


5.23 




120 


5.57 


5.66 


5.63 




NGC 3766 


3403.790 


3870.224 


4225.662 




2 
8 


3.24 


4.07 


4.11 
6.44 




16 


4.28 


4.26 


4.39 




23 


9.22 


9.37 


9.17 




25 


3.18 


5.30 


1.54 


• • • Be transient 


26 






5.46 




27 


5.38 


5.69 


5.59 




28 






6.38 




29 






4.47 




31 


3.99 


0.57 


— 6.42 


• • • Be transient 


33 






4.86 




36 


8.46 


8.56 


8.43 




40 






5.80 




41 


4.88 


5.28 


5.02 




42 


6.07 


6.02 


5.85 




45 


4.98 


5.10 


4.35 




46 






5.49 




47 


-11.33 


-11.29 


-6.43 


Be Be 


49 


4.37 


4.89 


4.70 




52 






4.61 




54 


5.18 


5.61 


5.31 




55 


5.03 


5.12 


5.25 




57 


5.67 


5.86 


5.66 




61 


4.09 


4.20 


4.24 


Be? 


69 


2.89 


3.19 


3.16 




72 


8.27 


8.65 


7.08 




73 


0.42 


3.28 


0.06 


Be? Be transient 


77 


4.61 


5.72 


5.50 




83 


4.07 


0.84 


2.93 


Be? Be transient 


87 
89 






5.35 
5.70 




92 


2.31 


2.73 


3.79 


Be? Be transient 


94 


5.94 


5.91 


5.86 


Be? 


96 


6.39 


6.87 


7.53 




98 


3.10 


0.21 


4.14 


Be transient 


101 


7.40 


7.44 






103 






5.21 




109 






5.10 




110 






4.83 




111 






4.56 




118 


7.40 


7.03 


6.73 




119 


-5.48 


-3.19 


-8.67 


• • • Be transient 


121 






6.54 




126 


7.46 


8.49 


7.74 


Be? 


127 


-4.95 


-4.00 


-5.03 


Be Be 


129 


8.66 


7.94 


7.96 




130 


2.00 


5.04 


5.09 


Be? Be transient 


131 






5.30 
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Table 2 — Continued 



Cluster, 
MG ID 


HJD (t - 2,450,000), 
WHa (A) (2005 Feb) 


HJD (t - 2,450,000), 
Who, (A) (2006 May) 


HJD (t - 2,450,000), 
Who, (A) (2007 May) 


Emission status Emission status 
in Paper 2 in this work 


133 


-1.15 


-9.14 


-12.72 


Be transient 


137 






6.02 




138 






5.09 




139 


5.46 


4.24 


5.05 


Be? Be transient 


142 






6.57 




145 






7.43 




148 






4.84 




149 






7.55 




150 






6.83 




154 


-40.07 


-33.56 


-33.40 


Be Be 


155 


7.92 


7.90 


7.42 




161 


3.92 


4.19 


4.18 




162 


9.15 


9.23 


8.85 




164 






6.46 




170 


4.81 


5.26 


5.34 




171 






5.11 




173 


7.79 


7.92 


7.80 




175 


6.67 


6.74 


6.64 




176 






3.54 




178 


6.32 


7.71 


7.40 




179 






4.25 




184 






5.72 




190 


6.81 


6.86 






192 






6.87 




196 


4.32 


4.39 


1.69 


• • • Be transient 


197 


7.65 


8.85 


8.49 




198 


-43.40 


-53.71 


-53.96 


Be Be 


200 


-6.66 


-4.62 


-5.08 


Be Be 



NGC 4755 3403.863 • • • 4223.717 

1 6:31 6:63 

6 9.50 • • • 10.33 

15 4.53 ••• 4.96 

18 7.42 • • • 2.57 • • • Be transient 

25 3.59 ••• 4.08 

27 9.94 ••• 11.11 

30 9.61 ••• 

31 4.85 • • • -2.09 • • • Be transient 
43 7.86 • • • • • • Be 

45 6.02 • • • 6.51 

47 5.48 • • • 6.42 

53 2.31 • • • 2.69 Be? 

54 9.35 • • • 10.06 
71 5.01 ••• 5.80 
73 8.33 ••• 9.11 
75 5.95 • • • 6.56 
93 5.85 • • • 6.79 

103 4.26 • • • 5.35 

120 5.20 • • • 5.99 

121 4.42 • • • 6.23 

129 -48.80 ••• -53.35 Be Be 

131 6.03 ••• •••Be? 

144 8.54 ••• 9.00 

148 6.48 ••• 6.99 

152 10.70 ••• 11.37 

166 10.81 • • • 11.60 Be? 

168 9.51 • • • • • • Be? 

169 4.99 • • • 
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Table 2 — Continued 



Cluster, 


HJD (t - 2,450,000), 


HJD (t - 2,450,000), 


HJD (t - 2,450,000), 


Emission status Emission status 


MG ID 


WHa (A) (2005 Feb) 


Whc (A) (2006 May) 


Whc (A) (2007 May) 


in Paper 2 in this work 


171 


— 55.57 




— 55.68 


Be Be 


172 


3.50 




4.20 


_' ' _' ' 


180 


— 44.23 




—36.48 


Be Be 


183 


4.11 




4.52 


_' ' _' ' 


187 


— 14.19 




— 13.86 


Be Be 


195 


2.47 




3.12 




205 


5.53 








215 


5.58 








222 


8.27 




8.23 


Be • • • 


223 


11.81 




11.08 


rse .' • • • 


252 


5.08 




' ' ' 




261 


8.13 




8.47 








o8o9.8z8 


4223. 81o 




2 




5.29 


4.78 


Be? Be transient 


4 




8.42 


8.14 




15 




6.44 


5.94 




32 




4.37 


4.38 




38 




8.38 


8.20 




45 




5.75 


5.53 




55 




5.62 


5.35 




57 




8.42 


7.77 




63 




7.15 


7.16 




67 




6.05 


5.59 


Be 


79 




— 5.46 




Be Be 


86 




8.24 


7.43 




90 




3.82 


3.17 


Be transient 


91 




5.18 


4.97 




92 




9.82 


8.32 




110 




4.57 


4.48 




111 




4.27 


4.22 




114 




6.48 


5.74 




116 




4.46 


4.49 




121 






/ .U4 




124 




4 49 


4 43 




132 




3.98 


4.03 




137 




4.12 


4.04 




156 




8.94 


8.38 




162 




3.14 


3.23 




165 




3.27 


3.39 




168 




3.94 


3.95 




169 




6.34 


5.99 




178 




6.87 


6.98 




185 




7.82 


7.12 




190 




3.52 


3.49 




213 




6.32 


6.19 




218 




4.14 


4.11 




221 




6.71 


6.96 


Be? Be transient 


235 




7.06 


6.97 




245 




5.65 


5.67 




268 




7.49 


7.96 




NGC 6664 




3870.933 


4223.930 




5 






11.20 




13 




6.88 


6.82 




33 




8.09 


7.84 




35 




6.55 


6.28 




45 




10.20 


10.32 
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Table 2 — Continued 



Cluster, 


HJD (t - 2,450,000), 


HJD 


(t - 2,450,000), 


HJD (t - 2,450,000), 


Emission status 


Emission status 


MG ID 


Whc (A) (2005 Feb) 




(A) (2006 May) 


Whc (A) (2007 May) 


in Paper 2 


in tliis work 


53 






8.19 


8.26 






62 






7.24 


7.12 






68 






8.90 


8.89 






73 






5.82 


5.79 






79 






3.82 


3.71 






81 






12.57 


12.17 






119 






6.42 


6.14 






126 






1.76 


2.84 




Be 


137 






-4.32 


-4.37 


Be? 


Be 


141 






11.41 


11.15 






143 






9.50 


9.00 






158 






11.35 


11.45 






168 






7.82 


7.59 






177 






-7.27 


-7.37 


Be 


Be 


182 






1.47 


1.74 


Be 


Be 


183 






7.20 


7.21 






207 






6.85 


6.73 






219 






5.69 


5.56 






221 






-18.49 


-18.68 


Be 


Be 


229 








10.48 






237 






7.13 


6.84 







^W^Ha may be contaminated by nebular emission in NGC 3293 and IC 2581. 



Table 3 
Fraction of Be Transients 



Cluster 


Number of Be 
Stars, n 


Number of Transient 
Be Stars, t 


Ratio of Transients 
to Total Be, p (%) 


CoUinder 272/Hog! 


y 16 6 


4 


67± 19 


IC 2581 


3 


1 


33 ±27 


NGC 3293 


5 


2 


40 ±22 


NGC 3766 


16 


11 


69± 12 


NGC 4755 


6 


2 


33± 19 


NGC 6231 


4 


3 


75 ±22 


NGC 6664 


5 





± 18 



Table 4 
Cluster Parameters 





Photometric fitting 


Eib-y) 


{V-Mv)o 


log age 




Cluster 


constant, C (mag) 


(mag) 


(mag) 


(yi-s) 


Reference 



CoUinder 272 
Hogg 16 
IC 2581 
NGC 3293 
NGC 3766 
NGC 4755 
NGC 6231 
NGC 6664 



0.49 
0.43 
0.51 
0.31 
0.27 
0.43 
0.44 
0.67 



± 0.06 
± 0.08 
± 0.05 
± 0.13 
± 0.02 
± 0.02 
± 0.06 
± 0.06 



0.34 
0.31 
0.31 
0.19 
0.17 
0.28 
0.35 
0.56 



11.55 
11.00 
12.00 
12.24 
11.6 
11.61 
11.07 
10.70 



7.11 
7.05 
7.14 
7.0 
7.16 
7.0 
6.9 
7.66 



Paper 2 
Paper 2 
Paper 2 

|Dufton et~aL\ j2006| 
Paper 'Z, Paper '6 



Dufton et aTlJloOe 



Perry et ar| l [l\jijl^ , Sana et alT] | |20Q6[ l 



"Paper"? 
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Fig. 1. — Ha profiles of Be stars, labeled by MG 
ID number. Spectra from 2005, where available, 
are shown with dotted lines, spectra from 2006 
with dashed lines, and spectra from 2007 with 
solid lines. 
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Fig. 2. — Ha profiles in the same format as Fig. 




Fig. 3. — Ha profiles in the same format as Fig. 
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Fig. 4. — The theoretical logarithmic relation- 
ship between Wna and y-Ra color {dotted line) 
is shown with the mean observed T4^Ha f^^nd the 
adjusted y-Ha color from Paper 2 for the tran- 
sient Be stars {filled diamonds), stable Be stars 
{open diamonds), and non-emission B-type stars 
{crosses). 



Fig. 6. — Cumulative Vsini distribution for each 
cluster, including both Be and normal B-type 
stars. 




400 



E{b-y) 



Fig. 5. — The fitting constants, C, are plotted 
against the clusters' reddening values, E{b — y). 
The best fit linear relationship is also shown {dot- 
ted line). 



Fig. 7. — Cumulative Vsini distribution for Be 
stars with stable disks and those with transient 
disks. The distribution for normal B-type stars is 
also shown for reference. 
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Fig. 8. — Tcff and logg are plotted for the cluster 
NGC 3293 with the evolutionary tracks of Schaller 
et al. (1992; solid lines). The ZAMS mass of each 
evolutionary track is labeled along the bottom. 
Normal B-type stars are shown as open squares, 
stable Be stars as open diamonds, and Be tran- 
sients as filled diamonds. 



Fig. 10. — Teff and logg are plotted for the cluster 
NGC 4755 in the same format as Figure [Sj 
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Fig. 11. — Toff and log g are plotted for the cluster 
Fig. 9.— Toff and logg are plotted for the cluster NGC 6231 in the same format as Figure § 

NGC 3766 in the same format as Figure [8] 
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